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ABSTRACT 


This  effort  is  intended  to  stimulate  further  research  and  development  of  an 
extremely  low  frequency  (ELF)  plasma  antenna.  The  study  explores  the  advantages  of  a 
transmitting  ELF  plasma  antenna  for  submarine  applications  and  draws  comparisons 
between  this  antenna  and  those  currently  in  use.  Three  promising  designs  are  discussed, 
and  mathematical  models  used  in  this  study  are  presented. 
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ELF  PLASMA  ANTENNA 


1.  INTRODUCTION 


The  extremely  low  frequency  (ELF)  plasma  antenna  concept  grew  primarily  from  five 
areas:  (1)  the  corona  mode  ELF  antenna,  the  properties  of  which  depend  on  ionization  of  the  air 
around  a  vertical  ELF  antenna  suspended  by  a  balloon  in  air;  (2)  the  Naval  Research  Laboratory 
(NRL)  reflector  plasma  antenna,  which  is  currently  being  tested;  (3)  a  1966  patent  entitled  “Laser 
Beam  Antenna”;  (4)  a  1965  report  (prepared  by  the  American  Nucleonics  Corporation  (ANC)) 
on  the  plasma  antenna  as  a  receiver  in  a  nuclear  environment;  and  (5)  mathematical  proof  of 
farfield  radiation  by  oscillating  plasmas. 

This  effort  draws  on  knowledge  and  experience  gained  in  these  five  areas  to  conceive  and 
design  a  transmitting  ELF  plasma  antenna  that,  if  implemented,  would  have  a  tremendous  impact 
on  submarine  ELF  communications:  the  convenience  of  antenna  transportability  would  replace 
unwieldy  vertical  electric  dipole  ELF  antennas,  balloons  suspended  in  air,  and  inefficient 
horizontal  electric  dipole  ELF  antennas. 

In  this  study  various  designs  were  explored  to  develop  an  ELF  plasma  antenna,  and  at 
least  three  designs  show  promise.  One  system,  which  would  use  the  reflective  properties  of 
plasmas  to  redirect  a  radar  signal,  is  under  development  by  NRL.  NRL  has  devised  a  plasma 
sheet  (which  currently  is  mechanically  rotated)  to  reflect  a  high-frequency  signal  radiated  by  a 
driving  antenna  (see  figure  1).  The  device  is  designed  primarily  for  surface  ship  applications.  A 
future  system  could  steer  the  plasma  sheet  electronically,  which  would  result  in  a  fast, 
multifunctional  antenna  reflector. 


Anode 


Figure  1.  NRL  Plasma  Reflector 
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Another  potential  application  of  the  reflective/transmission  properties  of  the  plasma  is  the 
reduction  of  the  radar  cross  section  of  an  antenna.  These  properties  can  reduce  the  size  of  the 
radar  cross  section  as  long  as  the  plasma  antenna  is  operating  at  a  frequency  below  that  of  search 
radars.  The  reduction  might  be  minimal  because  the  mounting  structure  usually  reflects  more 
radar  signals  than  the  actual  antenna  element,  but  even  some  reduction  is  important  for 
submarine  application  if  the  antenna  is  mounted  on  top  of  the  mast  or  combined  as  a  conformal 
antenna  on  a  stealth  sail. 

The  reduced  cross  section  is  probably  of  greater  importance  to  the  surface  ship 
community,  where  the  antennas  tend  to  be  relatively  large  and  can  contribute  significantly  to  the 
size  of  the  ship  cross  section.  For  example,  in  a  case  where  the  plasma  antenna  is  transmitting  at 
a  frequency  of  30  MHz  and  is  scanned  by  a  radar  operating  at  3  GHz,  the  amount  of  reflected 
energy  would  be  on  the  order  of  0.047  percent.  With  no  physical  structures  supporting  the  ELF 
plasma  antenna,  its  radar  cross  section  will  be  zero  because  its  natural  resonance  frequency  will 
always  be  less  than  the  radar  frequency;  the  antenna  can  therefore  be  operated  in  regions  where 
undetectability  is  desirable. 

Currently,  the  operating  ELF  antenna  is  the  horizontal  electric  dipole  (HED)  antenna. 
However,  horizontal  ELF  antennas  are  extremely  inefficient  and  must  be  located  where  large 
regions  of  low  ground  conductivity  exist.  The  vertical  electric  dipole  antenna,  with  and  without 
corona,  is  much  more  efficient  but  is  aerostat-supported,  unwieldy,  and  subject  to  “blowdown,” 
which  causes  this  tethered  antenna  to  assume  the  form  of  a  catenary.  The  vertical  electric  dipole 
ELF  antenna  is  experimental  and  has  not  been  implemented. 
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2.  ELF  PLASMA  ANTENNA  CONCEPT 


The  ELF  plasma  antenna  concept  uses  as  a  current  carrier  an  ionized  column  of  air, 
created  by  one  or  more  lasers,  with  a  length  on  the  order  of  the  vertical  electric  dipole  (with  or 
without  corona)  ELF  antenna.  The  vertical  length  can  be  12,500  ft,  as  it  was  for  the  vertical 
electric  dipole  ELF  antenna;  5.2  km,  as  in  the  corona  mode  antenna;  or  ideally,  as  in  some  of  the 
proposed  designs,  all  the  way  to  the  ionosphere  (30  km  to  70  km).  In  several  designs,  this 
ionized  column  of  air  is  forced  to  oscillate  at  ELFs. 

One  of  the  most  promising  designs  uses  a  high-powered  laser  to  ionize  a  pencil-thin 
column  of  air  and  to  control  the  oscillation  by  momentum  transfer  of  laser  photons  to  charge 
carriers  in  the  ionized  column,  which  affect  an  upward  current  by  momentum  transfer  and  a 
downward  current  by  laser  controlled  gravitational  relaxation.  An  alternate  method,  using  pure 
laser  control  of  the  level  of  ionization  and  current  oscillation,  is  to  use  electro-optic  modulation 
of  the  ionizing  laser  beam  to  send  an  ELF  wave  along  the  ionized  column.  The  electro-optic 
modulation  can  be  achieved  by  passing  the  laser  beam  through  an  electro-optic  crystal,  with  an 
ELF  oscillating  voltage  applied  transversely  to  the  crystal.  By  ionizing  a  thin  column  all  the  way 
to  the  ionosphere  and  sending  an  ELF  oscillation  along  this  column  by  electro-optic  modulation, 
this  oscillation  can  be  reflected  from  the  ionosphere  and  a  resonance  ELF  antenna  can  be  created. 

The  development  of  some  designs  has  examined  the  possibility  of  an  external  electric 
field  oscillating  at  ELFs  and  driving  the  plasma  column.  This  could  be  achieved  by  using  an 
oscillating  electric  field  from  an  antenna  at  the  base  and  focusing  the  electric  field  with  an 
electrostatic  lens.  Other  designs  have  included  the  use  of  acoustic  waves  to  drive  the  charge 
carriers  in  the  plasma. 


2.1  PLASMA  ANTENNA:  A  RECEIVING  DEVICE 

Research  performed  by  the  ANC  in  1965  (sponsored  by  the  Defense  Atomic  Support 
Agency)  demonstrated  the  ability  of  a  plasma  column  to  act  as  a  receiving  antenna  (ANC,  1965). 
The  objective  of  the  project  was  to  develop  an  antenna  to  measure  the  radiated  electric  field 
produced  by  a  nuclear  explosion.  Alternate  antenna  designs  were  sought  because  traditional 
metallic  antennas  were  adversely  affected  by  the  ionizing  radiation  produced  by  the  nuclear 
event.  The  antenna  prototype  developed  by  ANC  is  essentially  a  propane  torch  seeded  with 
cesium  or  barium.  The  additives  provide  the  flame  with  a  variable  conductive  component. 
Because  the  flame  is  conductive,  electromagnetic  (EM)  signals  couple  to  the  plasma  and  are 
received  via  a  cathode  follower  pickup  wire  located  within  the  flame.  The  wire  attracts  electrons 
that  give  rise  to  a  received  voltage.  The  received  wire  is  connected  to  a  traditional  spectrum 
analyzer  to  measure  the  induced  voltage. 

Measurements  performed  on  this  design  showed  that  the  propane  plasma  behaved  like  a 
whip  antenna.  The  frequency  response  was  linear  between  50  kHz  and  50  MHz.  The 
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researchers  suggested  that  the  low-frequency  response  of  the  antenna  could  be  improved  by 
increasing  the  conductivity  of  the  flame  (e.g.,  increasing  the  cesium  concentration). 

The  sensitivity  of  the  propane  plasma  antenna  was  measured  and  compared  to  the 
performance  of  a  14-in.  whip  antenna.  The  results  of  these  measurements  are  illustrated  in 
figure  2,  which  shows  the  received  signal  (terminal  voltage)  of  the  plasma  and  whip  antennas  as 
a  function  of  the  generator  voltage  at  the  driving  antenna;  the  plasma  antenna  appeared  to  be 
more  efficient  than  the  whip.  This,  however,  could  be  attributed  to  the  size  of  the  flame,  which 
might  have  been  larger  than  the  14-in.  whip  antenna,  resulting  in  a  larger  capture  area  (the  size  of 
the  flame  was  not  specified).  A  portion  of  the  study  was  also  devoted  to  developing  an 
equivalent  circuit  model  of  the  ANC  plasma  antenna;  the  final  model,  illustrated  in  figure  3, 
provides  a  relationship  between  the  resistance  and  capacitance  of  the  plasma  and  the  received 
voltage. 


Figure  2.  Comparison  of  Reception  of  ANC  Plasma  Antenna  vs.  14-in.  Whip  Antenna 
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Cp=  17.5  pF 


R2,  Cj  =  Burner  Assembly  to  Ground 


Figure  3.  Equivalent  Circuit  Model  of  Plasma  Antenna 


2.2  IONIZATION  LEVELS  AND  ENERGY  REQUIREMENTS  FOR  CREATION 
OF  AN  IONIZED  COLUMN  IN  THE  ATMOSPHERE  BY  LASERS 

The  primary  chemical  composition  of  air  is  shown  in  table  1 .  Since  the  density  of  air 
decreases  with  increasing  height,  the  corresponding  energy  required  to  ionize  the  column  will 
also  decrease  at  different  heights.  The  atmospheric  density  at  different  altitudes  is  listed  in 
table  2  (Roussel-Dupre  and  Miller,  1993b). 


Table  1.  Composition  of  Air 


Element 

Molecular  Weight 

Percent  by 
Volume 

Percent  by 
Weight 

First  Ionization 
Potential  (eV) 

Nitrogen 

N2  =  28.0 

78.09 

75.55 

14.53 

Oxygen 

O2  =  32.0 

20.95 

23.13 

8.5 

Argon 

Ar  =  39.944 

0.93 

1.27 

15.755 

Carbon  Dioxide 

C02  =  44.010 

0.03 

0.05 
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Table  2.  Atmospheric  Density  vs.  Altitude 


Height  (m) 

Density  (kg/m^) 

0 

1.225 

800 

1.156 

1600 

1.069 

2400 

0.967 

3000 

0.891 

4000 

0.819 

For  this  effort,  only  the  ionizing  of  a  volume  of  nitrogen  is  considered.  The  ionization 
potentials  in  eV  and  joules  are  illustrated  in  figures  4  and  5.  The  energies  shown  reflect  the 
values  required  to  reach  the  first  ionization  potential.  The  energy  required  to  achieve 
correspondingly  higher  ionization  states  can  be  computed  by  substituting  the  appropriate 
ionization  potentials  to  the  above  equations.  The  charts  can  be  used  in  conjunction  with  the 
parameters  of  existing  lasers  in  order  to  estimate  the  size  of  the  ionized  spheres  (e.g.,  when  the 
lasers  are  used  to  ionize  a  spot).  The  densities  in  table  2  were  substituted  into  the  previously 
described  equations  to  estimate  the  energy  required  to  ionize  a  1-cm  volume  of  nitrogen  at 
varying  heights.  The  results  are  illustrated  in  figures  6  and  7  for  both  eV  and  joules. 


Figure  4.  Ionization  Potential  (eV)  vs.  Volume  for  Nitrogen 
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Energy  (eVJ 


Figure  5.  Ionization  Potential  (J)  vs.  Volume  for  Nitrogen 


j 

Figure  6.  Ionization  Energy  (eV)  for  1  cm  vs.  Height  Above  Ground 
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Figure  7.  Ionization  Energy  (J)  for  1  cm3  vs.  Height  Above  Ground 


2.3  AIR  IONIZATION  OF  THE  PLASMA  COLUMN 


The  approach  for  computing  the  ionization  potential  is  to  determine  the  number  of  atoms 
in  a  given  volume.  The  ionization  potential  for  the  volume  is  then  computed  as  the  product  of 
the  number  of  atoms  times  the  ionization  potential,  which  can  be  done  using  the  equation  for  the 
speed  of  molecules  ( V)  (Roussel-Dupre,  1993a),  given  as 


where  p  -  pressure  (1 .01  x  105  N/m  =  1  atm)  and  r  =  density  (kg/m3 ).  This  can  be  arranged  to 
solve  for  the  density  r,  which  yields 


P  = 


3p 

F2rms  ’ 
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Velocities  for  nitrogen,  oxygen,  and  air  are  (Halliday  and  Resnick,  1986): 


Molecule 

Kmi  (m/s) 

n2 

493 

02 

461 

Air 

485 

The  computed  densities  for  nitrogen  and  oxygen  are  1.2466  kg/m3  and  1.426  kg/m3, 
respectively.  These  results  compare  favorably  with  the  tabular  values  provided  in  the  CRC 
Handbook  of  Engineering  (N  =  0.0013  g/cm3  and  O  =  0.0014  g/cm3).  The  number  of  atoms  per 
unit  volume  is  computed  as 


23 

where  NA  =  Avogadro’s  number  (6.023  x  10  )  and  As  =  atomic  weight. 

Converting  the  values  to  the  appropriate  units  yields  the  following: 

N2  =  2.68  x  1019  atm/cm3 
02  -  2.684  x  1019  atm/cm3. 

The  energy  required  to  ionize  the  volume  is  given  as  e(ev)  =  Vim  Nm,  where  Vi0B  is 
ionization  potential. 


9/(10  blank) 


3.  PRODUCTION  OF  ELF  CURRENTS 


The  particle  picture  of  laser  photons  shows  that  a  photon  has  momentum  and  can 
exchange  that  momentum  with  solid  matter.  Hence,  it  is  possible  to  drive  ions  or  electrons  in  a 
given  direction  by  colliding  them  with  photons.  This  is  the  idea  behind  creating  current  in  a 
plasma  with  the  use  of  a  laser  as  a  driver.  The  photons  can  be  made  to  be  absorbed  selectively  by 
the  ions  and  driven  through  inelastic  collisions  in  the  upward  direction.  The  laser  can  do  this 
over  one-half  the  ELF  period,  after  which  the  rapid  relaxation  of  the  ions  and  electrons  to  their 
respective  Maxwell-Boltzman  distributions  takes  place.  The  relaxation  time  period  exceeds  an 
ELF  period,  but  the  laser  can  be  used  to  slow  the  relaxation  time  down  to  half  an  ELF  period  by 
transferring  the  right  amount  of  momentum  from  the  photons  to  the  ions.  Hence  with  the  use  of 
the  laser  exclusively,  an  ELF  current  can  be  created  in  the  ionized  column.  (The  relaxation  time 
calculations  and  the  momentum  transfer  equations  are  provided  in  the  appendix.)  Figure  8  shows 
the  design  and  physical  concepts. 


Ionosphere 


Figure  8.  Laser-Induced  ELF  Plasma  Antenna  with  Current  Produced  by 
Photon  Momentum  Exchange  and  Maxwellian  Relaxation 
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3.1  LASER  AS  AN  IONIZATION  SOURCE 


Modulation  of  laser  beams  by  electro-optic  crystals  is  a  well-known  technology  that  can 
be  used  to  make  the  ionized  column  of  air  oscillate  at  ELFs.  In  this  method  a  laser  beam  passes 
through  an  electro-optic  crystal  with  electrodes  attached  to  a  voltage  source  that  oscillates  at  the 
rotational  frequency  (RF).  This  causes  the  laser  to  carry  an  ELF  wave  that  forces  the  charge 
carriers  in  the  plasma  to  oscillate  at  the  ELF  and,  hence,  to  radiate  an  ELF  signal.  To  create  a 
resonance  ELF  antenna,  this  wave  must  be  reflected  at  the  endpoint  of  the  plasma  antenna,  which 
occurs  because  the  natural  resonance  frequency  of  the  ionosphere  exceeds  the  ELF.  This,  in  turn, 
causes  the  plasma  ELF  wave  to  reflect  off  the  ionosphere  and  travel  to  Earth,  producing  a 
resonance  antenna.  Figures  9  and  10  illustrate  this  process. 


Plasma  Column 


Figure  9.  Use  of  the  Ionosphere  to  Give  Resonance  Standing  Wave  Plasma  Current 


12 


Beam 


Electro-optic 


Figure  10.  Laser-Produced  Resonance  Columns  on  the  Plasma  Column 

3.2  EXTERNAL  ACOUSTIC  AND  INTERNAL  ION-ACOUSTIC  DRIVEN 
PLASMA  COLUMN  ELF  CURRENTS 

Acoustic  waves  have  the  same  physical  properties  (for  the  same  frequencies)  as  EM 
waves,  but  they  have  much  shorter  wavelengths.  Acoustic  oscillations,  therefore,  are  able  to 
drive  the  charge  carriers  in  a  plasma  at  ELFs  over  a  shorter  distance  than  are  EM  oscillations. 
This  eliminates  the  need  for  an  aerostat-supported  ELF  transmitter,  but  the  ion-acoustic  antenna 
must  be  thicker  in  order  to  ensure  a  similar  dipole  moment  as  the  aerostat-supported  ELF 
antenna.  The  plots  in  figure  14  give  the  cross-sectional  area  of  these  types  of  antennas  if  the 
efficiency  of  the  antennas  is  100  percent — the  conversion  of  input  acoustic  energy  to  output  EM 
energy  is  done  without  losses  to,  say,  Bremstralung.  If  the  plasma  is  excited  externally  by  an 
acoustic  wave,  the  acoustic  wave  becomes  an  ion-acoustic  wave  in  the  plasma.  The  ion-acoustic 
wave  is  a  longitudinal  pressure  wave  in  which  the  ions  provide  the  inertia  and  the  electrons 
provide  the  restoring  force;  hence,  the  ion-acoustic  wave  is  an  ion  oscillation.  At  the  ion 
resonance  frequency,  the  ions  have  much  more  charge  density  than  do  electrons  oscillating  at  the 


electron  resonance  frequency.  Consequently,  the  ions  oscillating  at  resonance  and  set  to  equal 
ELF  can  provide  greater  charge  movement  and  a  greater  dipole  moment  than  the  electrons  can. 

The  other  attractive  feature  of  the  ion-acoustic  design,  shown  in  figure  1 1,  is  that  the 
plasma  can  be  contained  in  a  structure.  The  ionization  can  be  controlled  by  three  methods: 
electrodes  at  each  end,  microwave  heating,  and  laser  ionizing.  There  is  much  more  control  over 
the  ionization  in  this  column  than  there  is  in  other  designs  because  it  is  much  shorter  and  can  be 
contained  in  a  tube. 


Ion  Movement 


Acoustic  Source 


Figure  11.  Ion-Acoustic  Plasma  Antenna 


The  acoustically  driven  plasma  antenna  design  uses  acoustic  waves  as  the  transport 
mechanism  for  the  electrons  and  ions.  The  basic  theory  of  the  system  is  that  the  longitudinal 
acoustic  waves  propagating  through  a  fluid  (or  plasma)  medium  result  in  pressure  gradients  that 
in  turn  result  in  particle  movement.  In  this  case,  the  particles  are  ions  and  electrons,  and  the 
movement  will  give  rise  to  a  radiated  EM  wave. 

The  plasma  column  will  be  significantly  shorter  than  if  the  transport  mechanism  is 
attributed  to  an  electric  field,  so  it  is  conceivable  that  the  ionization  can  be  generated  using  either 
lasers  or  radio  frequency  (vibrational  and  rotational  excitation).  The  conceptual  view  of  the 
system  is  illustrated  in  figure  12. 
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Plasma  Column 


Ionization  by 
Lasers  or  Other 
Sources 


Acoustic  Source 


3.33  m 


Figure  12.  Acoustic  Plasma  Antenna 


The  acoustic  source  is  essentially  a  low-frequency  signal  generator/amplifier  combination 
that  modulates  a  transducer.  In  this  configuration,  the  acoustic  source  generates  a  wave  that 
“pushes”  and  “pulls”  the  ions  and  electrons,  thus  creating  movement.  An  example  of  the  particle 
movement  in  a  fluid  medium  is  illustrated  in  figure  13,  which  shows  that  particles  in  a  fluid 
(plasma  in  this  case)  cycle  through  compression,  displacement,  and  re-refraction  as  the  acoustic 
wave  propagates  through  the  medium.  As  the  amplitude  of  the  acoustic  wave  increases  (e.g.,  at 
point  p'),  the  plasma  will  start  to  compress  to  the  right  of  point  p'. 


Figure  13.  Acoustic  Wave  Interaction  with  Plasma 


The  compressed  plasma  increases  the  pressure  on  the  region  of  the  compressed  particles, 
which  in  turn  results  in  a  pressure  gradient  in  the  -s  direction,  causing  the  particles  to  accelerate 
in  the  +s  direction.  The  acceleration  increases  until  the  peak  pressure  point  arrives,  after  which 
the  compression  (hence,  pressure)  subsides  and  velocity  of  the  plasma  begins  to  decrease.  The 
velocity  slows  to  zero  as  the  acoustic  wave  begins  the  second  half  of  its  cycle.  As  the  negative 
portion  of  the  waveform  increases,  the  backward  velocity  of  the  particles  begins  to  increase,  and 
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it  continues  to  increase  until  the  negative  peak  is  reached.  The  velocity  then  slowly  returns  to 
zero  as  the  ambient  line  is  crossed  (Roussel-Dupre  and  Miller,  1993a). 

The  effect  of  the  acoustic  wave  on  the  ions  in  the  plasma  column  is  shown  in  figure  1 1 , 
which  illustrates  the  particle  compression  and  subsequent  movement  of  the  ions  as  the  acoustic 
wave  propagates  up  the  column.  The  acoustic  wave  produces  both  forward  and  backward 
movement. 


Sample  calculations  to  determine  the  height  of  the  plasma  tube  have  been  performed. 
The  relationship  between  velocity,  wavelength,  and  frequency  of  an  acoustic  wave  is  given  as 

Kfa=Va, 


where  Va  =  acoustic  velocity  (333  m/s),  Xa  =  acoustic  wavelength,  and  fa  =  acoustic  frequency 
(100  Hz).  Assuming  an  ELF  of  100  Hz,  the  wavelength  of  the  acoustic  wave  is  given  as 


333m  /  s 

Aa  = - = 

100  Hz 


3.33  m. 


This  implies  that  as  long  as  the  plasma  column  is  at  least  3.33  m  long,  it  should  be  capable  of 
radiating  an  EM  signal  at  the  acoustic  frequency.  The  result  is  a  wavelength  that  can  be  realized 
with  a  relatively  short  plasma  column. 

Setting  a  standing  acoustic  wave  on  the  plasma  column  causes  the  ions  to  move  at 
approximately  the  acoustic  frequency  and  the  electrons  to  move  at  another  frequency  (to  be 
determined).  The  electrons  can  move  at  a  different  velocity  because  of  their  smaller  masses. 

The  next  issue  of  concern  is  the  signal  generated  by  the  plasma  column.  Assuming  a 
sinusoidal  wave,  the  acoustic  pressure  is  expressed  as  (Roussel-Dupre  and  Miller,  1993a) 

P  =  Ppk  cos (oX  -kz +  (/>), 
the  acoustic  particle  velocity  (ions)  as 


z 

r  -  —  P  k  cos  (at  -kz +  </>), 
y  pc 


and  the  acoustic  intensity  as 


pc  ' 


The  ions  oscillate  and  cause  radiation  at  the  acoustic  frequency. 
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The  acoustic  power  is  given  as 


P 2 

P  =  Acoustic  Power  =  —  A , 

pc 

where  A  is  the  cross-sectional  area  of  the  plasma  column. 

As  a  result  of  conversion  of  energy,  the  relationship  between  the  input  power  and  output 
power  of  the  antenna  is  given  as 

P2 

P*=~A  =  Pou'+LOS’ 

pc 

where  Pout  is  the  output  power  and  losses  are  to  be  determined  (Bremstraulung,  etc.). 

Numerically,  the  estimates  of  the  output  power  are  computed  using  the  acoustic  power 
equation.  The  quantity  P2  is  computed  as 

p]=Plf'  10“  , 
and 

P2 

pm,=-±A> 

pc 

where  P2f  =20*10 ,p  =  1.3  kg/m3,  and c  =  333  m/s. 

Figure  14  illustrates  the  computer  output  power  for  acoustic  levels  ranging  from  80  dB  re 
20  pPa  (approximately  equal  to  the  noise  from  a  vacuum  cleaner)  to  140  dB  re  20  pPa  (jet  engine 
noise  level). 
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Output  Power  vs  Acoustic  Power 


vacuum  Acoustic  Power  (dB  re  20  uPa)  jet  Engine 

Cleaner 


Figure  14.  Output  Power  Density  vs.  Cross-Sectional  Area  and  Input  Power 

Each  of  the  curves  displayed  represents  a  different  cross-sectional  area.  The  cross 
sections  for  the  computations  varied  from  0.01  m2  to  1  m‘2.  The  graph  shows  the  radiated  power 
of  the  corona  mode  antenna  (0.05  W),  and  the  crossing  point  shows  that  the  acoustic  antenna  will 
be  capable  of  meeting  radiated  power  requirements. 


3.3  ELF  PLASMA  ANTENNA  DRIVEN  BY  FOCUSED  ELECTRIC  FIELD 

In  the  focused  electric  field  design  an  antenna  at  the  base  of  the  ionized  column  is  used  to 
produce  an  oscillating  electric  field  to  drive  the  current  carriers  in  the  plasma  at  ELF.  An 
electrostatic  lens  can  be  used  to  focus  the  antenna  electric  field,  which  otherwise  would  spread 
out  too  much  and  lose  its  effectiveness.  A  high-frequency  plasma  antenna,  contained  in  a  tube  as 
an  array  of  Hertzian  dipoles,  has  been  modeled.  (The  mathematical  model  and  corresponding 
computer  output  are  given  in  the  appendix.)  Each  dipole  is  represented  as  an  oscillating  electron- 
ion  pair — a  microscopic  Hertzian  dipole.  The  effect  of  these  atomic  Hertzian  dipoles  is  summed 
up  over  the  axis  of  the  antenna  to  yield  a  net  radiative  field.  The  Hertzian  dipole  design  was 
developed  for  a  high-frequency  antenna  oscillating  at  the  resonant  plasma  frequency  in  the 
megahertz  region  and  contained  in  a  tube.  The  Hertzian  dipole  radiative  mechanism  can  be 
extended  to  the  ELF  range  when  the  oscillating  ELF  electric  field  is  focused  by  an  electrostatic 
lens  into  the  ionized  column.  Figure  15  shows  this  design. 
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Plasma  Column 


Oscillating  Hertzian 
Dipoles  at  ELFs 


High-Powered 
C02  Laser 


Focused  Electric  Fields  as 


Figure  15.  ELF  Plasma  Antenna  Driven  by  Focused  Oscillating  Electric  Field 
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4.  HORIZONTAL  ELF  ANTENNA 


The  horizontal  ELF  antenna  design  is  not  a  high  priority  because  it  is  a  horizontal  electric 
dipole  antenna  concept  (see  figure  16).  Because  of  the  larger  mass  of  ions,  the  current  is 
primarily  an  ion  current.  (The  mathematical  model  for  this  current  appears  in  the  appendix.) 

The  appealing  aspect  of  the  horizontal  ELF  antenna  design  is  the  strong  control  that  exists  over 
the  plasma  current  by  an  oscillating  magnetic  field  at  ELFs.  This  is  usually  called  a  drift  current, 
and  the  direction  of  current  motion  is  perpendicular  to  gravity  and  the  magnetic  field.  As  the 
magnetic  field  changes  direction  at  ELFs,  the  current  oscillates  at  ELFs  with  a  large  dipole 
moment,  since  it  is  primarily  ion  current  oscillating  at  the  plasma  frequency  set  equal  to  the  ELF. 


Horizontal  ELF  Antenna 


Helmholtz  Coils 


I 


Gravity 


Ion  Current  Oscillating  at  ELF 


Figure  16.  Horizontal  ELF  Antenna 


The  drift  velocity  of  charged  particles  due  to  gravity  is  given  as 


-a  SxB 

VDG  ~  tO 


qa 
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where  a- e  for  electrons  and  a  — i  for  ions.  From  examination  of  this  equation,  the  following 
conclusions  can  be  drawn: 

•  Ions  and  electrons  drift  in  opposite  directions. 

mi 

•  Ions  have  a  velocity  that  is  —  >  velocity  of  electrons . 

Ttle 

•  Ions  drift  from  electrons,  causing  a  charge  separation  and  an  electric  field. 

Let  %  oscillate  at  frequency  w  =  100  kHz,  then 


B  =  Re  Beia* , 


B  ma  £xRe  Beja 

VDG  ~  C> 


B z 


V^G  =Re 


ma  gxReB 
d,  B2 


Therefore,  the  drift  velocity  oscillates  at  the  magnetic  field  frequency,  and  the  electrons  and  ions 
oscillate  in  opposite  directions,  providing  current. 
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5.  CONCLUSIONS 


All  the  designs  proposed  in  this  report  have  the  potential  for  transmitting  and  receiving 
radiation  patterns — even  though  the  physical  mechanisms  for  each  case  are,  to  varying  degrees, 
different.  The  ELF  plasma  antenna  models  and  the  high-frequency  Hertzian  dipole  model  are 
intertwined.  The  high-frequency  plasma  antenna  model  and  the  more  difficult  ELF  model  should 
be  developed  concurrently;  insights  gained  from  the  high-frequency  model  work  will  be  useful  in 
development  of  the  ELF  model. 

The  plasma  ELF  antenna  has  a  tremendous  advantage  over  the  horizontal  electric  dipole 
ELF  antenna,  which  lacks  portability  and  efficiency.  The  vertical  electric  dipole  ELF  antenna  is 
supported  by  a  balloon  at  least  12,000  feet  high,  is  extremely  cumbersome,  and  lacks  portability. 
Furthermore,  a  balloon-supported  metallic  antenna  is  not  stealthy.  The  ELF  plasma  antenna  is 
portable,  more  efficient  (because  of  its  extensive  use  of  focused  laser  light),  and  stealthy. 

This  report  is  intended  to  stimulate  research  and  development  of  the  ELF  plasma  antenna. 
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APPENDIX 


SUPPORTING  CALCULATIONS 


FARFIELD  RADIATION  OF  AN  OSCILLATING  PLASMA: 

ROUSSEL-DUPRE  MODEL 

The  nonlinear  interactions  occur  to  the  dynamics  of  the  plasma  itself.  An  analysis  by 
Roussel-Dupre  and  Miller  has  shown  that  the  normalized  radiated  electric  field  for  a  plasma 
cloud  (or  column)  is  given  as 


~  j  ,1  f  s  'j  ij 

E(s,x,z)  =  e~x /D  J  dz’—  e0a0(s,z')e  *'  +elal(s,z,)e 


A-l 


where  D  is  the  length  of  plasma,  and  s  is  the  frequency. 

The  magnitude  of  the  field  is  normalized  to  VqBq/c,  where  V0  is  the  initial  plasma 
velocity,  B0  is  the  background  magnetic  field  strength,  and  c  is  the  speed  of  light.  Therefore,  the 
magnitude  of  the  radiated  wave  is  directly  influenced  by  the  magnetic  field  and  the  plasma 
velocity. 

Roussel-Dupre’s  work  has  also  shown  that  the  radiated  EM  wave  will  have  the  same 
frequency  as  the  plasma  wave.  Further,  in  their  work  the  plasma  cloud  or  column  behaved  as  a 
dipole,  producing  the  pattern  illustrated  in  figure  A-l. 
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Figure  A-l.  Radiation  Pattern  of  Plasma  Column 

Numerical  analyses  performed  by  Roussel-Dupre  and  Miller  have  produced  the  values  for 
the  normalized  field.  Measurements  and  calculations  on  EM  emissions  from  plasmas  have 
varied  significantly.  Roussel-Dupre’ s  work  has  suggested  that  emissions  were  on  the  order  of  25 
mV/m  at  1.6  km.  At  the  other  extreme,  measurements  on  electrostatic  discharges  (1/2-in.  arc, 

26  A)  have  produced  radiated  electric  fields  on  the  order  of  150  V/m  at  1.5  m. 


RADIATED  FIELD  FROM  PLASMA  OSCILLATIONS:  MHD  MODEL 


Others  have  shown  mathematically  that  plasma  oscillations  are  also  responsible  for 
radiating  EM  fields.  Montgomery  and  Tidman  (1964)  have  shown  that  the  basic  cold  plasma 
equations  are 


di  d 
ck  dx 


(nv)  =  0 


3 


d  ^  i  m 
—  xE  =  - — — 
dx  c  <k 


—  •E  =  4m(N0-n)y 


d  1  dE 

rx5  =  — — 
dx  c  ck 


4m 

- - nv 

c 


d_ 

dx 


•  5  =  0, 


where  N0  is  the  ion  density,  n(x,t)  is  the  density  of  zero  temperature  electron  gas,  and  v(x,t)  is  the 
velocity.  Montgomery  and  Tidman  (1964)  use  a  perturbation  scheme  to  develop  a  set  of  coupled 
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equations.  The  equations  are  solved  to  determine  the  radiated  field,  and  the  solution  is  divided 
into  first-  and  second-order  solutions.  The  first-order  electric  field  is  (Ginsburg,  1964): 

d<f> 

E  =  — sin (aoet)  +  AT  sin (vt-K^x  +  a). 

The  second  term  of  this  equation  represents  the  transverse  component  of  the  electric  field. 
The  velocity  of  the  wave  is  given  as 

v2  =K2c2  +a>', 

where  K  •  At  =  0.  The  corresponding  magnetic  field  is  given  as 
c 

B  =  —  KxAt  sm(vt- K»x  +  a) . 
v 

The  source  size  of  the  plasma  oscillation  must  also  be  defined  and  is  given  as 
#  =  4>0(l  +  S*x)e~xl,l} , 


where  d  is  the  deviation  from  spherical  symmetry,  and  L  is  the  thickness  of  the  plasma  slab. 
The  magnetic  field  at  a  distance  x  from  the  plasma  is  then  given  as  (Ginsburg, 1964) 

32V2 mco)'  x 

where  n  is  the  gas  density  function  of  x  position  and  time,  e  is  the  electron  charge,  m  is  the 
electron  mass,  Nq  is  the  ion  density,  (f>  is  the  source  size,  coe  is  the  electron  frequency,  coe  is 


* 4nN0e2S\ 


V  m  ) 


,  and  k  is  the  wave  number  = 


CD. 


Since  the  electric  field  is  related  to  the  magnetic  field  via  E  =  rjH,  it  is  obvious  that  the 
radiated  electric  field  is  also  a  function  of  the  plasma  electron  frequency.  The  frequency  of  the 
second-order  magnetic  field  is  equal  to  twice  the  plasma  frequency  and  is  valid  when  coeL  <c. 

Similar  equations  can  be  obtained  using  the  plasma  ion  frequency  instead  of  the  electron 
frequency.  The  ion  frequency  is  believed  to  have  a  more  pronounced  effect  because  the  ion  mass 
is  significantly  greater  than  an  electron  mass. 
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ELF  PLASMA  ANTENNA  PRELIMINARY  CURRENT  PREDICTIONS 


A  1-kW  carbon  dioxide  laser  with  a  beamwidth  of  2  mm  produces  3.18  billion  w/m2  of 
energy  intensity.  If  this  converts  to  pure  radiation  intensity  in  the  farfield,  then  the  energy 
intensity  in  the  laser  is  equal  to  the  farfield  Poynting  vector.  The  farfield  Poynting  vector  is 
equal  to  the  square  of  the  H  field  times  the  characteristic  impedance.  The  characteristic 
impedance  is  numerically  377.  Setting  the  square  of  the  H  field  times  the  characteristic 
impedance  equal  to  the  laser  intensity,  we  can  solve  for  the  H  field  and  get  56,000  A/m.  Using 
Amperes  law  and  choosing  somewhat  arbitrarily  the  ELF  distance  as  4000  m  as  the  length  of  the 
antenna,  we  set  the  product  of  the  H  field  and  the  antenna  length  equal  to  the  current  carried  by 
the  plasma  antenna.  We  obtain  several  hundred  amps,  which  is  a  far  greater  current  than  can  be 
carried  by  a  metal  antenna.  Of  course,  this  result  depends  on  100-percent  conversion  of  the  laser 
energy  into  the  farfield  ELF  plasma  radiation;  however,  if  the  energy  conversion  is  only  a  few 
percent,  this  is  much  more  current  than  can  be  generated  on  a  metallic  antenna. 


RELAXATION  TIMES 

When  the  atmosphere  is  ionized  with  lasers,  the  electrons  and  ions  relax  to  separate 
Maxwellian  distributions  because  of  the  difference  in  kinetic  energy  and  gravitational  potential 
energy.  The  Maxwellian  distributions  for  electrons  and  ions  (Roussel-Dupre,  1993b)  are  given 
as: 


Electrons: 


e-{m.y*+m.g A)/2W 


Ions: 


Ji  \2nkTJ 


3/2 


-(m(vJ+m,g/i)/2*r 


The  relaxation  times  for  electrons  and  ions  (Roussel-Dupre  and  Miller,  1993b)  are  given 


as: 


Electrons: 


m/\2kTe ) 


x. « 


3/2 


8;zwee4[d>(l)  -4e~l  /  V^r]  In  A 
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Ions: 


_ m^jlkTf2 

%mteA  [o(l)  -  4e~'  /  Vtt]  In  A 

where 


and 


_  1/2 
Jc3T3  1 

tzw  y  ze2  ’ 

The  net  current  density  from  relaxation  time  is 


From  Krall  and  Trivelpiece  (1973,  p.  306),  the  relaxation  times  for  gas  discharges,  with 
densities  on  the  order  of  1014  cm'3  (similar  to  those  computed  for  our  ionized  plasma  column), 

n  <7 

are  on  the  order  of  2  x  10' s  for  electrons  and  2  x  10' s  for  ions.  This  suggests  that  the  transient 
current  generated  is  on  the  order  of 


1014  1014 
J~  2»10-9  +  2»  10-7  * 

The  result  is  obtained  by  multiplying  the  above  equation  by  1.6  x  10'19  coulombs.  The 
current  per-unit  volume  is  then  8,080  A/cm3.  The  current  density  is  then  computed  by 
multiplying  the  volume  by  a  length.  Assuming  a  1-cm  length  (i.e.,  the  length  between  spheres), 
the  current  density  is  8,080  A/cm  .  The  results  of  this  analysis  suggest  that  each  ionized  sphere 
must  be  revisited  on  the  order  of  10’9  s. 


The  pulse  durations  for  some  common  lasers,  in  the  mode  locked  operation,  are  as 
follows: 


Laser  Type  Pulse  Duration  (s) 

He-Ne  6x  10'10 

Nd:YAG  7.6  xlO'11 

Ruby  1.2  x  10'11 
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When  ionization  occurs,  te  =  lO^s  is  the  relaxation  time  for  electrons  to  the  Maxwell- 
Boltzman  statistics  of  electrons.  During  this  time  the  electrons  move  in  the  z  direction  with 
respect  to  the  air  molecules,  coming  from 

<*><  -<*>*, 


where 


\zpeen'z,kT' 

^\=±x - , 

jpte~mrZ/kTr 

0 

and 

z 

\zpabfrm~zncr~ 

\pairem'*z,kT~ 

0 

where  re  is  the  density  of  electrons  at  sea  level,  and  rair  is  the  density  of  air  at  sea  level.  The 
relaxation  time  for  ions  is  te  @  10"7s  and  will  not  be  included  in  this  approximation. 

me  =  9.1  lxlO'31  kg,  mass  of  an  electron, 
r=300K, 

Te  =  electron  temperature  (to  be  determined),  and 
mair  -  the  weighted  mass  of  an  air  molecule. 

Hence  the  electron  current  that  exists  over  time  te  @  10'9  s  is 

If  a  pulsed  laser  can  be  synchronized  with  the  recombination  time  te  such  that 

2| 

x  ,  =  —  «  X  +  X  . 
ciV  e  c 

the  mode  number  N  and  laser  length  /  to  give  “continuous”  current  can  be  determined. 
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ZERO  RADAR  CROSS  SECTION  OF  ELF  PLASMA  ANTENNA  COLUMN 


One  of  the  added  benefits  of  the  plasma  antenna  is  a  reduced  radar  cross  section.  This  is 
important  to  submarines  if  the  antenna  is  mounted  on  top  of  the  mast  or  combined  as  a  conformal 
antenna  on  a  stealth  sail,  but  it  is  probably  of  greater  importance  to  the  surface  ship  community, 
where  the  antennas  tend  to  be  relatively  large  and  can  contribute  significantly  to  the  size  of  the 
ship  cross  section. 

This  section  describes  the  methods  by  which  EM  waves  are  transmitted  and  reflected  by 
the  plasma  antenna.  The  generic  antenna  configuration  is  illustrated  in  figure  A-2. 


z=  o 
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A 

Incident  Wave 


p’  Reflected  Wave 

1  : 


Plasma  Region 


Dielectric  Constant 
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t* 

y 


Transmitted  Wave 


Figure  A-2.  Wave  Reflection  at  Vacuum/Plasma  Interface 


The  boundary  between  the  vacuum  and  the  plasma  region  is  depicted  at  Z  =  0.  The 
electric  and  magnetic  field  components  of  the  incident  wave  are  given  as 
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r 


E,=E Br=E,0e'**, 

where  EM  is  the  amplitude  of  the  incident  wave.  The  components  of  the  reflected  wave  are 
given  as 

E'X  =  REX  0e+ik°z  By  =  REx0e+ik°z , 

where  R  is  the  reflection  coefficient.  The  transmitted  wave  components  are 


E,  =  TEe 


+ikoZ 


Eke 

xop 

CO 


where  T  is  the  transmission  coefficient. 


The  relationships  between  the  wave  numbers,  frequencies  in  a  vacuum  and  in  the  plasma, 
are  given  as 

kgC2  =  co 2  (vacuum) , 


k2c2  =  co2 


1- 


co, 


2\ 


0)2  J 


(plasma) , 


where  kp  is  the  wave  number  in  plasma,  ko  is  the  wave  number  in  a  vacuum,  wp  is  the  plasma 
frequency,  and  w0  is  the  frequency  in  vacuum. 

The  reflection  coefficient  is  given  as 

K  -k„ 

D  _9 _ P_ 

k0  +  k/ 

and  the  transmission  coefficient  as 


*o  +  V 

Using  the  relationships  between  the  plasmas  and  the  vacuum,  the  conditions  listed  in 
table  A-l  are  evident.  Table  A-l  shows  that  for  high-frequency  waves,  there  will  be  perfect 
transmission.  Thus,  as  long  as  the  plasma  antenna  is  operating  at  frequencies  below  that  of 
search  radar,  the  radar  cross  section  of  the  plasma  antenna  will  be  less  than  an  equivalent  metallic 
antenna.  As  an  example,  consider  the  case  where  the  plasma  antenna  is  transmitting  at  a 
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frequency  of  30  MHz  and  is  scanned  by  a  radar  operating  at  3  GHz.  The  amount  of  reflected 
energy  would  be  on  the  order  of  0.047  percent. 


Table  A-l.  Reflection  and  Transmission  Coefficients  vs.  Frequency 


Frequency 

R 

T 

Comment 

High  (w  »wp) 

0 

1 

Perfect  transmission 

w*wp 

0<i?<l 

o<r<2 

Partial  transmission,  partial 
reflection 

w  =  wp 

1 

2 

Oscillations 

Low  (w  «wp) 

-1 

0 

Perfect  reflection 

HERTZIAN  DIPOLE  ARRAY  MODEL  OF  PLASMA  ANTENNA 
WITH  COMPUTER  OUTPUT 

The  radiated  electric  field  produced  by  a  line  of  charges  in  the  plasma  is  calculated  by 
modeling  electron/ion  pairs  as  Hertzian  dipoles.  The  total  radiated  field  is  then  the  summation  of 
the  fields  radiated  by  each  individual  dipole. 


The  force  on  an  electron  in  a  time-varying,  harmonic  electric  field 


(3* 


is  given  as 


r  r 
F  =  ~eE  ’ 

where  e  =  1.6x10'  C.  The  force  is  also  expressed  as 

r  d2x2  2  r 
F  =  m~ff  =  m(0  x’ 

where  m  is  mass,  and  co  is  the  angular  frequency. 

The  dipole  moment,  N& p,  is  equal  to  the  charge  times  the  distance  between  the  charged 
particles,  as  illustrated  in  figure  A-3. 


©  G 

Figure  A-3.  Distance  Between  Charged  Particles 
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Mathematically,  the  dipole  moment  is  given  as  =  qrf.  The  dipole  moment  per  unit 
volume  p  is  given  as 


where  N  is  the  number  of  electrons  in  the  plasma  per  unit  volume: 


r  _  Ne2  r 
P  =  ~1mo2E' 

D  is  given  as 

r  r  r  r  Ne2  r 

d=£oe+p=£oe-^e- 


Simplifying  the  equation  yields 


r 

D-£o 


1-- 


O) 


CO 


r 

E ’ 


where  cop  =  ^  frequency. 

In  this  application,  the  Hertzian  dipoles  are  formed  by  two  charged  particles,  as  illustrated 
on  the  left-hand  side  of  figure  A-4.  At  a  time  t  =  774,  the  current  flows  between  the  charges  until 
the  charge  on  the  particles  is  reversed  (t  =  Til).  The  current  then  changes  direction  as  the 
charges  again  reverse  themselves  between  the  two  particles. 


Hertzian  Dipole 


©  O  ©  O 


i  i 
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f=774  f=772  <=3774 


Figure  A-4.  Hertzian  Dipole 
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The  insertion  loss  (IL)  product  for  these  miniature  dipoles  is  given  as 


Az  =  jcap, 


where  p  =  qDz.  The  electric  and  magnetic  field  components  for  a  Hertzian  dipole  are  given  as 


klAze 


rjkr 


4n 


-sin#  , 


and 


r 


H~ 


,  klAze~Jkr  .  „ 


The  electric  and  magnetic  fields  are  perpendicular  to  each  other.  The  wave  impedance  is  given 
as 


The  length  between  the  charged  particles  is  computed  by 


900  MHz  =  ^ 
2tc 


n(1.6»10“19)2 


9.1 1  •  10~31)(8-85»  10~12)  ’ 


where  n  =  10'8  electrons/m3.  The  linear  spacing  of  the  electrons  is  then  given  as  spacing  =  yfn  . 
Substituting  values  yields  an  approximate  spacing  of  1  mm.  Thus,  Dz  =  1  mm. 

If  the  antenna  is  1  cm  long,  and  operating  at  900  MHz,  then  the  sum  of  the  Dz  from  -0.5 
cm  to  0.5  cm  should  yield  the  desired  pattern  for  the  single  line  of  ion/electron  pairs.  Based  on 
the  coordinate  system  found  in  figure  A-5,  the  field  is  given  as 

sin0|  , 


Pi  V  s  4 nr 
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r 


and  the  Poynting  vector  is  then  given  as 


to 


r  r 
EXH 


j]  ^&|/|AzV 
2  v  Anr  J 


sin2  6  . 


Figure  AS.  Coordinate  System  for  Plasma  Model 
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